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ABSTRACT 

We investigate a thermo-resistive instability in the outer crusts of 
magnetars wherein a perturbation in temperature increases ohmic 
heating. We show that magnetars of characteristic age Tagc ^ 10"' yr 
are unstable over timescales as short as days if strong current sheets 
are present in the outer crust. This instability could play an impor- 
tant role in the thermal and magnetic field evolution of magnetars, 
and may be related to bursting activity in magnetars. 



1 INTRODUCTION 

Soft gamma repeaters (SGRs) are a class of highly magnetized neutron stars (magnetars) that exhibit 
persistent X-ray emission, interrupted periodically by short bursts of gamma rays. The more frequently 
occurring short bursts have typical durations of 0.01 - Is, with peak luminosities of ~ 10**^ erg s^^. Giant 
flares are much more rare and energetic outbursts, with peak luminosites 10^ — 10^ times larger than 
the short bursts. Thre e giant flares have b een observed to date: from SGR 05 26-66 on March 5, 1979 



( Helfand fc LongHl979[ ) 



Mazets et al 



i Wld)). SGR 190 0+14 on August 27, 1998 jHurlev et al.ll 19991 ). and 



SGR 1806-20 on December 27. 2004 jflurlev et al.lboosl '). SGR 1806-20 has produced the largest observed 



giant flare, releasing ~ 10*® erg of energy over 380 s. Precursors to two of the three giant flares have 
been identified. Th e 2004 giant flare was preceded by a Is long event, 140 s prior to the initial hard spik e 



(jHurlev et al 



20oi). A similar precursor to the August 27 giant flare was observed ([Ibrahim et al.ll200ll ). 
with duration ~ .05 s preceding the initial spike by 0.4 s. 

SGR flares ar e thought to represent the r elease of magnetic energy, though the trigger mechanism 



Thompson fc DuncanI jl995h propose that the short bursts observed from SGRs are the 



remains uncertain. 

result of fracturing the crust of the neutron star by magnetic stresses. Recen t calculations indicate that 
neutron star crusts fail catastrophically under stress ([Horowitz fc Kadaull2009l b The rigid crust could act 
as a gate, releasing magnetic energy when magnetic stresses cause the crust to fail. Larger events such as 
giant flares could be caused by large scale readjustment of the magnetic field due to magnetic instability 
([Thompson fc DuncanI Il995l ). As the field readjusts, magnetic reconnection lead s to the form ation of a 
pair p lasma which is injected into the magnetosphere ( Thompson fc Duncan|[T995 b Alternatively, Lvutikov 
([20061 1 argued that the short rise time of the 2004 giant flare requires that the magnetic energy be stored in 
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the magnetosphere rather than the stellar interior. In this scenario, slow untwisting of the internal magnetic 
field eventually leads to sudden relaxation of the twist in the magnetospheric field, releasing the energy 
necessary to power the fiare. 

The initial configuration and subsequent evolution of magnetic fields in highly-magnetized neutron 
stars is a complicated problem. The field evolves continu ously due to the efi'ects of ohmic decay, ambipolar 
diffusion, and Hall drift. Recentlv. lPons fc GeppertI l|2007l ) studied the evolution of magnetic fields in neutron 
star crusts, emphasizing the importance of Hall drift. Their results indicate that Hall drift of crustal fields 
can create small-scale magne tic field structures, and that those structures can drift to regions of higher 
resistivity. The simulations of iPons fc GeppertI (|2007l ) were restricted to magnetic fields in the inner crust. 

In this paper, we focus on the outer crust, and show that large currents can lead to a thermo-resistive 
instability, affecting the thermal evolution of the star. As the instability evolves, large portions of the crust 
may melt, allowing the magnetic field to evolve on a timescale much faster than the average ohmic decay 
and Hall timescales. The enhanced magnetic evolution resulting from instability may be related to flare 
activity in magnetars. 

This paper is organized as follows. In section 2, we describe the relevant physics of the thermo-resistive 
instability in neutron star crusts. In section 3, we describe the physical processes that lead to evolution 
of the neutron star magnetic field. In section 4, we describe our neutron star model and give details of 
calculations of the instability growth rate. Section 5 contains discussion and our conclusions. 



2 THERMO-RESISTIVE INSTABILITY 



The dipole fields of magnetars are inferr ed to be in t he ra nge ^ (0.5 — 20) x 10^^ G, based on spindown 
measurements of SGRs and AXPs ( see iMereghettil (|2008l ) for a review). In order for t he magnetic field 
to be stable in neutron stars, it must contain both toroidal and poloidal components (jMarkev fc Tavler 
19731 ). The toroidal component in a stable configuration has a twisted torus shape, and may be an order of 
magnitude larger than the poloidal component ( Braithwaite , 2009; ') . 

Large crustal currents associated with the toroidal field would produce significant ohmic heating due 
to the relatively high resistivity in the outer crust. Ohmic heating can account for the observ ed trend of 



surface temperature that increases with surface field observed in neutron stars with i? > 10 G ( Pons et al 



20071 ) . Cooling simulations show that a heating layer in t he outer crust, as wou ld arise from current decay. 



can explain the high surface temperatures of magnetars (jKaminker et al 



2006 



Current decay is determined primarily by electron-phonon interactions for a temperature T below the 
melting temperature Tmeit. In this regime, the electrical resistivity scales as T, so that a small increase 
in temperature leads to increased heat dissipation (Fig. 1). The additional heating raises the temperature 
further, and a temperature runaway may develop if thermal transport is unable to quench the feedback 
process. As we show in this paper, this instability can occur in the outer crusts of neutron stars, where the 
electrical resistivity is relatively high and the thermal conductivity is low. 



© 0000 RAS, MNRAS 000, 000-000 



Thermo -Resistive Instability in Magnetar Crusts 3 



s 



SO 
o 




7 8 9 10 

Log Temperature 

Figure 1. Resistivity of the crust at a density of 10^ g cm~^. The melting temperature at this density is ~ 3 X 10* K. 



3 OHMIC DECAY AND HALL DRIFT 

The evolution of the magnetic field in neutron stars after birth is determined primarily by ohmic decay and 
Hall drift. The ohmic decay timescale is Tohm ~ rj~^L^ , where L is the typical magnetic field length scale 
and 77 is the electrical resistivity. A typical value for the outer crust at temperature T — 10* K is 
L n2 / 



Tohm = 3 X 10 



Vlkm/ Vl 



n 



,10-3cm2s-i' 

For a magnetar of age IC* yr, crustal currents from the initial field should still be present. [We note that 
the Tohm was smaller early in the star's life, since the temperature and resisitivity were hi gher. 



Hall drift creates small scale magnetic structures in the crust over the Hall timescale (jPons fc Geppert 

2003), given by 



THall 



cB 



(2) 



where is the electron density and B is the magnetic field strength. A typical value for the outer crust is 



"^""~KtI^) (toSg) 



(3) 



Hall drift can concentrate currents in the crust. The induction equation, neglecting ohmic dissipation, is 



= -V X (^^-^—{V X B)x B 



(4) 



dt \4:Tvnee 

To illustrate how Hall drift may affect the magnetic field, consider a magnetic field in cylindrical coordinates, 
with only an azimuthal component which depends on r. 



B = B^{r)4>. 
Eq. (4) becomes 

9B_^yf ' B, d{rB,) ^ 
dt Wvrnee r dr 



(5) 



(6) 



Since the magnetic field has only r dependence, and the quantity inside the parenthesis is in the f direction, 
the curl is zero, and the field is stationary. 
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Outward drift of the field can replenish currents in the outer crust that are decaying through ohmic 
diffusion. In order to g et an outward drift of the field, the field must have z-dependence, as considered by 
( Pons fc GeppertI 120071 ). This would correspond to a field strength that varies from the magnetic pole to 
the magnetic equator. The induction equation for the (/>-component of the field gives 



dB^ 2 dB^ 



(7) 



dt THaii r dz 

If dB/dz is positive, the field in the crust will increase over a typical timescale of THaii- 

We assume henceforth that strong currents exist in the crust at an age of ~ 10** years, and explore the 
consequences. 



4 CALCULATIONS 

4.1 Equations and Boundary Conditions 

We model the outer crust as an infinite slab, with x pointing into the star (see Fig. 2). The thermal evolution 
of the neutron star outer crust is described by the energy conservation equation, 

c4p, T) ^ = i^!ZklI) j2 ^ ^ . [^(^^ (^)j ^ 

where p is the density, c„ is the specific heat, j is the electric current, is the neutrino emissivity, and k 
is the thermal conductivity. The magnetic field evolution is described by the induction equation 

^^ = -Vx (r,(p,r)VxB(x)), (9) 
where the magnetic field is related to the current by 

= X B{x). (10) 

As justified below, we work in an approximation in which magnetic induction can be ignored, so we need 
not specify boundary conditions on B. The slab consists of of 3 regions - an atmosphere with no magnetic 
dissipation, an outer crust, and an isothermal inner crust/core. The atmosphere extends from a density of 
10'*gcm~^ at the stellar surface to lO^gcm"^, t he outer boundary of the crust. Fo r the atmosphere and 
outer crust zones, we employ the density model of lFriedman fc Pandharipandd (|l981f ). The inner crust /core 
zone is assumed to be an infinite heat reservoir, beginning at a density of p = lO^^gcm"'^. 

Boundary conditions. At the stellar surface we choose the unperturbed temperature to be, 

T^Ts, (11) 

and require that the heat flux at the boundary equal the blackbody emission rate at temperature Ts, 

.f^.T/. (12) 

As a simple model of the toroidal component of the neutron star magnetic field, we introduce to the 
outer crust zone a current sheet of width L: 

j(:r)=joe-'"~^"/'^''5, (13) 

where jo is the amplitude of the current, and a^o the location of the peak current. This analytic form 
allows a large, nearly uniform current near the heating peak, falling off sharply for \x — xo\ > L (Fig. [H]). 
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Figure 2. Neutron Star model. The shading represents the region of ohmic heating 
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Figure 3. Sample current sheet and associated magnetic field. The current peaks at xq = 100 m, approximately at 
the center of the outer crust. 



The magnetic field resulting from the current lies in the y-z plane, with variation in the x direction. In our 
models, the heating region is near the center of the outer crust, with characteristic width much smaller than 
the crust thickness to ensure negligible heating at the boundaries. We note that this current model leads 
to large pressu re gradients i n the crust. To form a stable current model, a complex geometry is required, 



such as that of 



thermo-resistive instability. 



Braithwaitd (120091 ). We use the simplified current model described here to demonstrate the 
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Figure 4. The specific heat at 10 K. The solid line is ionic specific heat, and the dashed fine is electronic specific 
heat. 



4.2 Input Physics 



For th e electrical and thermal conductivities in the crust, we use the analytical expressions of iPotekhin 
(|l999f ) for electron-ion collision frequencies in Fe matter. We ignore the effects of the magnetic field on the 
conductivities. Small jumps in the conductivities occur at T = Tmeit, which we smooth to avoid numerical 
problems. The heat capacity of the crust has contribut ions from ions an d relativistic electrons. The specific 
heat due to ions for solid matter (7 > 150) is given by Ivan Ripeii (|l99ir ). 



where 

22.75^2 



1 + 



log (7) 



7 ■ 



log (150) 



1/3 



erg cm K , 



(14) 



(15) 



4 = 5.4 X 10' 



Tg erg cm K ^ , 



(16) 



is the Coulomb plasma parameter, pa is the density in units of IQ^gcm"^, Tq is the temperature in units 
of 10^ K, Z is the ionic charge, and A the atomic weight. The contribution from relativistic, degenerate 
electrons to the specific heat is 

■.no J 

where no = 0.16 fm^'^ is the nuclear saturation density. The specific heat of ions and electrons are shown 
in Fig. 4 at 10* K. 

In strong ly-magnetized neutro n star crusts, the neutrino lumino sity is dominated by neut rino synchro- 
ton emission (jAguilera et al.ll2008l ). The rate of emission is given bv lBezchastnov et al.l (|1997l ) 

= lO'^Bi^rl ergcm"^s^\ (17) 

where _Bi5 is the magnetic field strength in units of 10^^ G, and Tg is the temperature in units of 10^ K. 
The ratio of neutrino emission to ohmic heating H is given by 



= 7 X 10"^ B?5 tI 



{ ^ V 

V 10-* cms/ 



10^^ (erg cm--' s~ 



2-11/2 



(18) 
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Figure 5. A sample temperature profile, with Ta = 10^ K, Bmax = 5 X 10^^ G. The dashed curve is the melting 
temperature of the lattice. Shading indicates the heated region. This model indicates that a portion of the crust is 
molten, but the heated region is solid. 



For the range of temperatures and magnetic fields we consider, the energy lost through neutrino emission 
is negligible compared to the ohmic heating rate and we neglect it in this analysis. 

Using these boundary conditions, we integrate eq. ([U to the boundary of the outer crust. An example 
of a typical equilibrium crust temperature profile is plotted in Fig. [5] corresponding to the current model 
shown in Fig. |3] The heat current for a typical crust model is shown in Fig. |6l Most of the energy flux is 
into the star, in the direction of increasing thermal conductivity. The energy is then lost to neutrinos from 
the core. 



4.3 Stability Analysis 



We now examine the stability of the equilibrium state . We perform a stability analysis of the outer crust 
using eq ((Sjl, substituting T{t) — To + STe''^ , where 5T is the perturbation mode and 7 is the growth (decay) 
rate. For the range of heating models we consider, instability growth rates are fast compared to magnetic 
diffusion and we neglect evolution of the magnetic field in our calculations. Appendix A contains further 
justification of this approximation. Neglecting magnetic evolution and neutrino emission, the perturbed 
energy balance equation is given by 



47r , .2 
Ct.7 2"^ Jo 



dTo 
dx dx 



dx"^ 



5T~ - 



dlio ^ ^ j dTp 
dx dx 



dST 
dx ' 



(19) 



where primes indicate differentiation with respect to T at fixed x. 

Perturbation mode boundary conditions. To determine the perturbation mode gradient at the outer 
boundary, (p = lO'' g cm~'^), we integrate through the atmosphere for several values of Ts to determine 
the dependence of the temperature gradient on the temperature. Because there is no ohmic heating in the 
atmosphere section, the temperature at the outer boundary Tot and the temperature gradient dTob/dx are 
uniquely defined for a given surface temperature. Therefore, dTob/dx is a function of Tob- 
dTob 



dx 



f{Tob 



(20) 
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Figure 6. Absolute value of the heat current for a sample crust model. The heat current flows towards the surface 
and the core from the heating peak. Most of the flux is lost to the core. 



Allowing perturbations to the temperature at the outer boundary for eq. (20) gives 
d{Tob + STob' 



= f{Tob + 5Tob). (21) 
tion f{Tob) is well behaved, eq (21) to first order in STob is 

^ f{Tob) + f{Tob)5Tob, (22) 



dx 

where primes indicate differentiation with respect to Tob- Using eq. (20) we can solve for f'{Tob), 
Keeping only perturbed terms of eq. (22), we arrive at the outer boundary condition. 



dSTob _ _d_ r dToi 
dx dTob \ dx 



STob. (24) 



We assume the inner crust/core of the neutron star to be an infinite heat reservoir. Therefore, at the 
inner boundary we require that the temperature perturbation vanish, 

ST,b = 0. (25) 



We solve eq. (|19|l simultaneously for the temperature perturbation mode 5T(x) and the mode growth rate 
7 for a given crust temperature profile and current distribution, given by eq. (9). We satisfy the set of 
mixed boundary conditions through a shooting algorithm. A sample temperature perturbation mode for 
the heating model presented previously is plotted in Fig. 7. The instability growth rate 7 as a function of 
the maximum field in the crust due to crustal currents is plotted in Fig. 8. 

To see the scalings of 7 with the parameters of the problem, it is useful to perform a local plane wave 
analysis to obtain an approximate analytical expression. Solving eq. (|19|) . assuming that k has only weak 
dependence on T and x and keeping only the dominant terms (determined by numerical experiment), we 
obtain 



d'^ST If 47r , .2 



dx'^ n 



Cvl j-q Jo 



ST. (26) 
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Figure 7. Sample unstable temperature perturbation mode, corresponding to the same crust model as Fig. 6. Dashed 
curve indicates the (normalized) electric current amplitude. Units are arbitrary. 
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Figure 8. Instability growth rate 7 vs maximum field in the crust, for heating peaked at p = 3 X 10^ gem . The 
minimum field required for instability is iJ = 3 X 10^^ G. 



Substituting plane-wave solutions for 5T with fc = 1/L gives an approximate expression for the mode growth 
rate, 



1 

7 ~ — 



47r , .2 Kp 

C2 '^"^^ L2 



(27) 



where all parameters are evaluated at the heating peak. Eq. (|27|l shows the competition between the ohmic 
heating (the first term), and cooling through thermal diffusion (the second term). Changing the surface 
temperature affects the growth rate by changing the crust temperature at the heating location. If the crust 
temperature is much greater than TmM , the temperature sensitivity of the resistivity rj becomes negligible 
and the heating feedback effect is lost, thereby stabilizing the system. Increasing the current amplitude 
gives more thermal energy to drive the instability. The growth rate is somewhat insensitive to the choice of 
L. However, the relationship between the magnetic field and the current (eq. (12)) indicates that for a fixed 
current amplitude jo, larger values of L correspond to larger magnetic fields. To obtain crust models with 
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realistic magnetic field amplitudes, the electric current must be concentrated in a relatively small region 
of the crust. Finally, the instability growth rate is highly dependent on the heating location xo due to the 
spatial variation of both the thermal conductivity and the electrical resistivity. A fixed current will produce 
the most heat in regions of high resistivity. Heat produced in regions of low thermal conductivity is most 
likely to lead to instability, as the heat is not efficiently carried away. For these reasons, heating near the 
surface of the star is most unstable. 

We calculate the instability growth rate for a wide range of heating models - characteristic heating 
widths range from 5 m < L < 50 m, and heating peak locations 10^ g cm~^ < po < 10^" g cm~^ . In this 
range, the current amplitude is virtually zero at the inner and outer boundaries. We consider current 
amplitudes corresponding to maximum magnetic fields from lO'^^ G to 10^^ G. The steady state surface 
temperature Ts is not the same for each heating model. We calculate crust temperature profiles with 10® 
K < Ts < lO'^ K, seeking models for which the temperature in the heating region is less than or equal to 
the melting temperature, and the temperature at the inner boundary is 10^ K < Tit < 10* K. Maintaining 
a crust temperature below Tmeit ensures that the resistivity feedback is operating at the heating location, 
and 10^ K < Tit < 10* K ensures that the core temp erature is in the appropriate range for a magnetar 



with a characteristic age of 10* yr ( Aguilera et al. 2003). These conditions lead to a narrow range of values 



for the surface temperature. Fig. 9 illustrates the region of instability for given values of the maximum 
field, Bmax, and the core temperature of the neutron star for fixed heating width and location. Subject 
to the conditions mentioned above, our results show that a minimum crustal field of S = 3 x lO'^^ G is 
required to destabilize the crust. Instability growth rates for several models are plotted in Fig. 10. To find an 
approximate expression for the critical field required to give instability (7 > 0), we substitute j ~ cBo/'i-nL 
in eq (|27|l and solve for the critical field, 

Bcot = 2.5 X 10^^ (ax lOiBergs'^-icm-iR-i) (iTTFfc^) ^'^^^ 
This estimate agrees well with our numerical solutions. 



5 DISCUSSION AND CONCLUSIONS 



Large crustal currents associated with the magnetic field of magnetars may lead to a thermo-resistive 
instability in the crust. Calculations of the instability growth time for a model neutron star crust give 
typical growth times of weeks to months. These timescales are short compared to the ohmic diffusion 
timescale of the magnetic field. 

We conclude that the instability identified in this paper may operate in neutron star crusts for a 
wide range of physical parameters relevant to magnetars. Heating may be located anywhere in the outer 
crust, while magnetic length scales need only be comparable to the crust thickness or smaller. Instability 
occurs for crust temperatures Tcruat ~ 5 x 10* K or lower, characteristic of relatively young magnetars. 



Tago ~ 10* yrs ( Aguilera et al.ll200a ). Th is result coincides with the inferred age of magnetar candidates 



associated with supernova remnants (see Mereghetti ( 2008h for a review). Additionally, we find that only 
heating models that produce large magnetic fields {B > 10^^ G) will produce instability, so this instability 
is specific to magnetars. The characteristic temperatures and magnetic fields at which the thermo-resistive 
instability occurs suggest and intriguing connection between this instability and magnetars. We note that 
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Figure 9. Diamonds represent values of the crust magnetic field and the neutron star core temperature for which 
unstable modes are found. The solid line indicates the approximate boundary between unstable and stable parameter 
space. All models used a heating width of 20 m, with heating location po = 3 X lO^gcm"''. 




Figure 10. The instability growth rate as a function of the neutron star core temperature and the crust magnetic 
field. All models used a heating width of 20 m, with heating location po = 3 X lO^gcm"'' . 



our simplified treatment of the current sheet is likely to overestimate the critical field required for instability 
by a factor of order unity. A stable current sheet necessarily has a more complicated structure than assumed 
here. The components of the current that we have neglected will lead to further heating. 

While we restrict our models to heating in the outer crust, instability in the inner crust could arise in 
the same way. Heat and charge transport mechanisms are no different, and we expect the scaling of eq. (|27p 
for the growth rate to apply to inner crust heating. However, because of the larger thermal conductivity in 
the inner crust, deeper crustal currents would have to be larger to produce similar instability growth rates 
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Model 


B (G) 


Q(ergs 1) 


Eb (erg) 


1 


1015 


2 X 1035 


5.7 X 10*5 


2 


5 X lOi"' 


1036 


1.4 X 10-*^ 


3 


10"' 


2 X 1036 


5.6 X lO'*'^ 



Table 1. Magnetic field, heat deposition rate and magnetic energy in the crust for 3 heating models. 

to those calculated here. As the heating is moved to deeper layers, the minimum magnetic field required 
for instability becomes greater than 10^^' G. 

The magnetic energy in the crust for our higher heating models is sufficient to power even the largest 
giant flares. The magnetic energy contained in the magnetic field in the crust is given by 

EB = 4n J r^^^dr, (29) 

where we integrate from the inner boundary to the stellar surface. We calculate the maximum field in 
the crust, the rate of energy deposition due to crustal currents, and the magnetic energy in the outer 
crust for several heating models (Table 1). The energy released during the largest giant flare to date was 
approximately 10*® erg. 

Future work is required to determine the nonlinear evolution of the magnetic field once an instability is 
triggered. Solving the coupled energy conservation and magnetic induction equations as a function of time 
will give insight into this problem. Our solutions for the steady state temperature indicate that a portion 
of the crust is molten. As the instability grows, higher density regions of the crust may melt, reducing 
the maximum magnetic stress that could be supported by the crust. Simulations of the magneto-thermal 
evolution could provide a link between the instability we have identified and glitch behavior in magnetars. 



Appendix A - Magnetic Induction 

In this Appendix we perform a local, plane-wave analysis of eqs. (19) and (20). The perturbed energy 
conservation equation including magnetic induction is given by 

dKo ^ ,dTo] dST \8n . 1 ^. 
ox ox ] ox icr J 

The magnetic induction equation (eq. (2)) can be written in terms of the current j, 

1^ = -V X (V X no)- (31) 
Upon linearizing the induction equation, we have 
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rio 



dx^ 



Sj - 

dx dx 



dx dx dx^ dx^ 



ST 



dST , . d^ST 
1130- 



dx 



dx^ 



(32) 



Magnetic induction can be neglected if the induction term is small compared to the ohmic heating 
term in eq. (31), integrated over the heating region: 



47r , .2 



— Vojofix), 



(33) 



where f{x) = Sj/ST. We find an approximate expression for f{z) using plane wave solutions for ST and Sj 
in eq. (33) and solving for Sj: 

-1 r 



Sj 



,2 , d r]o 



2ik 



^dri' djo , ,d'^jo , d^ri' . drf . ,djo ,. ,2 



ST. (34) 



Wc take for the wavenumber fc = 1/L, the characteristic width of the heating region. We approximate 
the induction term by evaluating the variables at the heating peaJj, 



8n 



Vojofix) 



877 ^ . 



(35) 



where L is the heating length scale, rjp is the resistivity at the heating peak, jmax is the maximum current, 
and f{x) is evaluated at the heating peak in the plane wave approximation. We approximate the heating 
term in the same way, 



i 



47r , ,2 ^ 47r , ,2 

30 — "^^hVpJmax-: 



(36) 



(37) 



where 77^ is ^ evaluated at the peak. Using (18) and (19), eq. (16) becomes 

i^hflp'Jmax ^ i^h f^p Jmax J [X ) ■ 

For the range of magnetic field models we consider, the heating term is much larger than the induction 
term. We plot the ratio of the heating term to induction term in Fig. 11 for several models. We see that 
magnetic induction is negligible for the large fields of interest. 

As a second test, we evaluate the ohmic decay time of the magnetic field by dimensional analysis. 
Dimensionally, the ohmic decay time is 

where rj is the resistivity at the heating peak, and L is a characteristic length scale. Ohmic diffusion 
can be neglected if the ohmic decay timescale is much longer than the instability growth timescale, Td ^ Tg. 
Fig. 12 shows the ratio Td/Tg for several heating models. 
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Figure 12. Ratio of instability growth timescalc to ohmic decay timescale for several heating models. Magnetic 
induction can be neglected for Tgjr^ <C 1. 
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